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We present the results of the first principle calculations of the energy bands in graphene and 
their symmetry classification. The valence bands and four lowest conduction bands are classified 
according to their symmetry at the points T and K. Merging of the bands is interpreted in the 
framework of the group theory. 
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INTRODUCTION 

Since graphene was first isolated experimentally [lj, it 
is in the focus of attention of both theorists and exper- 
imentalists. Obviously, our understanding of graphene 
starts from the knowledge of its energy bands. It is com- 
monknowledge that the highest valence (n) band and the 
lowest lying conduction (jr*) band merge at the Fermi 
level at point K (the corner of the Brillouin zone). The 
dispersion law of these two bands (obtained in the tight- 
binding approximation) is also well known since the sem- 
inal paper by Wallace The other bands in graphene 
were also studied previously, although in a much less de- 
tailed way than the above mentioned two. In particular, 
there were performed first principle calculations of the 
band structure of graphene [4|-l8| . 

However, a crucial element in our understanding of the 
band structure of graphene is still missing, that is the 
symmetry analysis of the bands. Such analysis (for any 
periodic crystal Hamiltonian) is most easily performed 
with the help of elementary group theory [§, [l(| ■ The 
symmetry of crystal is characterized by a point group R. 
Any operation of this group (save the unit transforma- 
tion) takes a general wavevector k into a distinct one. 
However, for some special choices of k some of the op- 
erations of the crystal symmetry group will take k into 
itself rather than into a distinct wavevector. These par- 
ticular operations are called the group of k; it is a sub- 
group of the full symmetry group of the crystal. Lines in 
the Brillouin zone for which the group of the wavevector 
contains elements other than the unit element are called 
symmetry lines. At special points in the Brillouin zone 
the group of the wavevector may be larger than that on 
symmetry lines which thread it; these are called symme- 
try points. We may use a state (states) corresponding 
to such a special wavevector to generate a representation 
for the group of k [i[l0|. (For an arbitrary wavevector, 
of course, the group of the wavevector is simply E, and 
the only irreducible representation which may be realized 
by the states corresponding to such a wavevector is the 
unit representation.) 



We may classify states at a wavevector corresponding 
to a symmetry point according to the irreducible repre- 
sentation of the group of the wavevector at that point. 
As the wavevector then moves away from the point, the 
group of the wavevector becomes smaller and some of 
the degeneracies will be split. We may determine the ir- 
reducible representations into which the original repre- 
sentation will split. Conditions relating the irreducible 
representations of adjoining points and lines are called 
compatibility relations and were discussed for the first 
time in Ref. [Il[. The early classical papers in the field 
are reprinted in Ref. [ijj]- The energy band calculations 
are most commonly carried out along symmetry lines in 
the Brillouin zone. 



SYMMETRY CLASSIFICATION OF THE 
ENERGY BANDS 

We have calculated the band-structure of graphene 
along the T-K-M-T line with the code Elk [l3( which im- 
plements the full potential linearized plane wave method 
(FP-LAPW) [3] and the local-density approximation 
(LDA) exchange-correlation (xc) potential. The k-vector 
grid was unshifted with 64 x 64 x 1 points with all the 
symmetries applied. The separation d between the peri- 
odically repeated in the perpendicular direction graphene 
layers was chosen at 200 bohr. With this, the elimination 
of the inter-layer interaction has been ensured. 

The results of the calculation are presented in Fig. Q] 
The curves marked have been found well converged with 
respect to the further increase of d and they represent 
the intrinsic nine lowest-lying valence and conduction 
bands of graphene. While the overall agreement with pre- 
vious pseudopotential-based band-structure calculations 
0-0, EH is found, one important point warrants a special 
noting: In Ref. [l5|, the two low-lying (at T point) bands 
were attributed to image-potential states. Since, in con- 
trast to Ref. IT5I we do not impose an additional image 
potential, and since LDA for xc potential does not ac- 
count for the latter [l6| , we conclude that the two bands 
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in question are ordinary Kohn-Sham LDA energy bands 
rather than the image-state bands. 

The gray background in Fig. [T] represents the contin- 
uous spectrum which is extrapolated in the d — > oo limit 
from the quasi-continuous one obtained from the calcu- 
lations. The edge of the continuum is at approximately 
4.6 eV at the T point which gives the vacuum level in this 
calculation. The bands entering continuum (red lines in- 
side the gray background) turn into resonances, i.e., the 
corresponding wave- functions leak out from the graphcne 
plane into vacuum. However, this leakage is apparently 
weak as indicated by the stability of the corresponding 
bands with respect to the d increase. The studies of the 
width of the resonances, including those corresponding 
to the higher bands, will be submitted for publication 
separately. 




Figure 1: (Color online) Graphene band-structure 
evaluated with use of the FP-LAPW method. The 
red-marked lines are well-converged single graphene 
layer bands, while gray background corresponds to 
continuous spectrum. 



The symmetry points in graphene are point T, the cen- 
ter of the Brillouin zone, points K, which are corners of 
the Brillouin zone, and points M , which are the centers 
of the edges of the Brillouin zone. The symmetry lines 
are T-K , T-M and K-M lines. The small group of vector 
k at point T is D§h- The small group of vector k at point 
K is D^h- The small group of vector k at lines T-K is 

The representation of the group D§h we can obtain 
noticing that 



D 6h = D e xd 



(1) 



Thus representation, say, A\ of group D$ begets two rep- 
resentations of the group D§h'- A\ g and A\ u , where the 
letter g or u means that the representation is even or odd 
with respect to inversion respectively The z-axis is al- 
ways chosen perpendicular to the graphene plane. When 
considering the T-K line, the x axis is chosen in the di- 
rection of the line. 
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Table I: Character tables 
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Table II: Character tables; e = e 27 ™/ 3 



The compatibility relations of the groups are obtained 
by taking into account that irreducible representations 
of the groups Dgh and D^h become reducible when we 
consider them as the representations of the group Ci v . 
The expansion of these representations with respect to 
irreducible representations of the group C^ v can be easily 
done using the tables of characters (Tables U and Hi]) . 
Equation 



9^ 



(2) 



shows how many times a given irreducible representation 
is contained in a reducible one. All the relevant compat- 
ibility relations are presented in Table IIIII 

Now let us come directly to the classification of the en- 
ergy bands. We assume that the representation realized 
by the lowest valence band at point T is the maximum 
symmetry representation A\ g . From Table IrTTl we imme- 
diately come to the conclusion that the band at the line 
r — K is characterized by the representation Ai . The rep- 
resentation Ai is compatible with one-dimensional rep- 
resentation A x and two-dimensional representation E' at 
the point K. But because of degeneracy in this point 
we come to the conclusion that the lowest valence band 
realizes the representation E 1 at point K, and the other 
band, merging with it at point K, at line T-K realizes 
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Table III: Compatibility relations for the honeycomb 
lattice (T-K line) 
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the representation B\. Hence the third band, merging 
with the second one at the point K, at line T-K realizes 
the representation A±. The second and the third band 
realize at point T cither the representation E\ u or E 2 g. 
Assuming that these two bands produce a bonding (that 
is symmetric) orbital, we come to the conclusion that 
they realize at the point T representation E\ u . Taking 
into account that Uh — IC2, we deduce that these valence 
bands are even with respect to the ah operation. 

Let us come to two bands which merge at the Fermi 
level at the point K . It is common knowledge that 
they are antisymmetric with respect to the Uh opera- 
tion (see Sec. J. Among the two bands in question the 
valence band (7r-band) is characterized at point T by the 
representation A\ u (bonding and, hence, symmetric or- 
bital), and the conduction band (7r*-band) is character- 
ized at point r by the representation Big (antibonding 
and, hence, antisymmetric orbital). At line T-K they 
are characterized by the representations A 2 and B2 re- 
spectively, and at point K they are characterized by the 
representation E" (HElII]. 11 

was mentioned already in 
the classical paper [17j | , that a two dimensional represen- 
tation realized by the wave functions at point K means 
that the dispersion law in the vicinity of this point can 
be presented by a circular cone. Similar reasoning can 
be applied to the lowest lying conduction bands. All the 
results of symmetry classification are presented on Fig. 

m 

The most prominent features in the band structure pre- 
sented in Fig. 1 is the merging of the bands. In the next 
Section we will interpret merging of the 7r and n* bands 
at point K and merging of the a bands at point T using 
the tight binding model. In the Section after the next 
one we will interpret merging of the bands at point K 
using the quasi-free electrons model. 

TIGHT-BINDING HAMILTONIAN + GROUP 
THEORY 

It is interesting to look at the proposed symmetry clas- 
sification of the bands in general and at their merging in 
particular from the point of view of the tight-binding 
model. In isolated form, carbon has six electrons in the 
orbital configuration ls 2 2s 2 2p 2 . When arranged in the 
honeycomb crystal, two electrons remain in the core Is 
orbital, and are traditionally ignored in band calcula- 
tions. The remaining four electrons occupy four valence 
bands; there are also conduction bands, of which a few 
lowest ones are of interest. The tight-binding Hamil- 
tonian for sp-bonded systems includes four orbitals per 
atom: s,p x ,p y ,p z . The Hamiltonian being symmetric 
with respect to reflection in the graphene plane, the 
bands built from the p z orbitals decouple from those 
built from s,p x ,p y orbitals. The former are odd with 
respect to reflection in the graphene plane, the latter are 



even. The classification with respect to reflection in the 
graphene plane already done, we can talk about Cq v and 
C% v groups instead of D§h and D^h- 

The n bands 

The structure of graphene can be seen as a triangular 
lattice with a basis of two atoms per unit cell, displaced 
from each other by any one (fixed) vector connecting two 
sites of different sub-lattices, say 5 = —a (1, 0). The gen- 
eral Hamiltonian for the n bands is 

TT _ _ ( Ea *'(a)e* k - a Ea *(a + S)e^+^ \ 
V£a**(a + <5K ik(a+5) Ea^'(a)e ika )' 

(3) 

where a is an arbitrary lattice vector, that is a linear 
combination of ai = | (3, y/3) , a 2 = § (3, — . 
The selection rule for matrix elements gives 

^(a + <5)e ,;K - (a+,5) = 0, (4) 

a 

where K is a corner of the Brillouin zone. In fact, we 
are dealing with the product of two functions: i(a + 5) 
realizes the unit representation of the point symmetry 
group C3 (the full symmetry group of the inter-sublatticc 
hopping is C^ v , but the restricted symmetry C3 will be 
enough to prove the cancellation). As far as the func- 
tion e lK '( a+l5 ) is concerned, rotation of the lattice by the 
angle 2n/3, say anticlockwise, is equivalent to rotation 
of the vector K in the opposite direction, that is to sub- 
stitution of the three equivalent corners of the Brillouin 

zone: Kx ->■ K 2 -> K 3 -+ K 1; where K : = (§f , ^) , 

K * = (°. - iSs) and K ^ = Thus due to 

the rotation e lK '( sl+s > is multiplied by the factor e 2 and 
realizes x — iy representation of the group C3. Because 
each of the multipliers in Eq. ((4]) realizes different irre- 
ducible representation of the symmetry group, the ma- 
trix element is equal to zero. Simply speaking, at a point 
K the sublattices become decoupled, and this explains 
the degeneracy of the electron states in this point (these 
points) or, in other words, merging of the two branches 
of the single Brillouin zone. 
On the other hand, generally 

£t'(a) e * Ka ^ 0. (5) 

a 

To understand this statement consider the maximum 
symmetry group of the intra-lattice hopping: C'ev The 
function i'(a) realizes the A\ representation of the group. 
Applying Eq. |(2J) we see that reducible representation of 
the group C$ v , realized by the two functions e lK a and 
e jK a can be decomposed as A\ + B-2- 
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In addition, the tight-binding model provides us with a 
simple explanation of why the dispersion law in the vicin- 
ity of the merging points is linear, that is why these points 
are Dirac points. The dispersion law for the Hamiltonian 
([5]) is given by equation 



F(E,k) = 0, 



(6) 



where 
F(E,k) = 



E a *(a + <5)e' ik -( a + 5 > 



£ a **(a + <5)e-' ik -( a + 5 ) E + J2 a t'(a.)t 



(7) 



In mathematics the Dirac points, we are dealing with, 
are called conical points of the surface; if the surface is 
given by Eq. ©, the conditions for the conical points 
are 



dF 
dE 



0. 



dF 



0. 



(8) 



Recalling the rule for differentiating of a determinant, 
we realize that Eq. (j4j guaranties that the conditions ((8]) 
for k = K. This explains linearity of the spectrum in the 
vicinity of the points K(K'). 



The a bands 

We will concentrate upon the dispersion law of the 
a bands along the T-K line. Mirror reflection in this 
line exchanges the graphene sublattices, changes y to — y, 
and does not change the wave vector. It means that 
the valence bands labeled by A\ are constructed from 
sub-lattice symmetric combinations of s and p x orbitals 
and from sub-lattice antisymmetric combinations of p y 
orbitals. The valence band labeled by B\ is constructed 
from sub-lattice antisymmetric combinations of s and p x 
orbitals and from sub-lattice symmetric combinations of 
p y orbitals. 

We can be more specific speaking about the valence 
states at points T and K . The states with the symmetry 
A\ g and A\ are the sub-lattice symmetric combinations 
of s orbitals. The states with the symmetry E\ u and E' 
are the (degenerate) sub-lattice symmetric combinations 
of p x and p y orbitals. 

Looking at the conduction bands we realize that the 
two bands which merge at the point T with the symmetry 
Eig are probably of the same nature as the two valence 
bands which merge at the point T with the symmetry 
Ei u , and the states with the symmetry E^ g are the (de- 
generate) sub-lattice antisymmetric combinations of p x 
and p y orbitals. (In the language of quantum chemistry 
sub-lattice symmetric and antisymmetric combinations 
are bonding and antibonding orbitals respectively.) 
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Table IV: Characters of natural representation 



QUASI-FREE ELECTRONS + GROUP THEORY 

If we consider separately a and ir states we can ignore 
the third dimension and treat graphene as purely two 
dimensional. The Hamiltonian of the electrons is 



H = Hn 



V. 



(9) 



where Hq is the Hamiltonian of the free electrons and 
V is the crystalline potential. The following reasoning 
uses the extended Brillouin zone scheme [22[ and is the 
application of group theory to quantum mechanics [20| . 

All the eigenstates of the Hamiltonian H with the 
same |k| are degenerate. The potential V having lower 
symmetry partially (nearly completely) removes the de- 
generacy. To find the maximum splitting we should find 
among the continuum of the vectors with the given value 
of |k| those which realize an irreducible representation of 
this or that sub-group of the space group of the crystal. 
For arbitrary value of |k| each state is symmetric only 
with respect to the sub-group of translations and real- 
izes a one-dimensional representation of this sub-group. 

However, the states with the wavevector K + b, where 
b is an arbitrary vector of the inverse lattice (corners of 
the Brillouin zones) have higher symmetry. They realize 
a representation (infinite dimensional) of the group which 
includes in addition to translations the point group C% v . 
We will expand this representation into the irreducible 
ones in two stages. 

First we will divide all the vectors K + b into the sets 
such, that elements of a given set arc connected with 
each other by the operations of the group C^ v . Such 
sets turn our to be triplets. For example, the states 
with the wavevectors Ki, Ki — bi, Ki — bi b 2 , where 
b x = |f (1, V3) and b 2 = (l, are the elemen- 

tary reciprocal-lattice vectors. Another example is the 
states having wavevectors Ki— b2, Ki+b2, Ki — 2bi b 2 . 
Still another example is two triplets having the same 
energy: K x — 2hx — 2b, K x + b^K! — b x + b 2 and 
K 1 -b 1 -2b 2) K 1 -2b 1 ,K l + b l +b 2 . 

At the second stage, applying Eq. |(3J) we decompose 
the reducible representation R of the symmetry group 
Cs v , generated by a triplet into the irreducible represen- 
tations. The triplet generates natural representation of 
the group C$ v [21|. Applying Eq. |2J) we obtain 



R = A 1 +E. 



(10) 



It should be noted that the result does not depend on V 
being in any sense small [20( ] . 

Recalling the symmetry of the states relative to reflec- 
tion in the plane of graphene, Eq. (|10[) can be rewritten 
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R = A[ + E' 



for a states, and as 



E" 



(11) 



(12) 



for 7r states. 

A picture that emerges becomes clear in the extended 
Brillouin zone scheme, where the Brillouin zones fill the 
whole plane. Each Brillouin zone being a hexagon, three 
zones meet at their corners; two of them merge. 

To find the dispersion law in the vicinity of the merging 
point we may consider a reduced Hamiltonian, 



H 



Hn H12 
H21 H22 



(13) 



where each element is a linear function of k x , k y (k is the 
deviation of the wavevector from the merging point), and 
we have shifted the energy axis, such that the energy in 
the merging point is now zero. The dispersion law for 
the Hamiltonian (|13|) is given by equation 



E — c%ki zb d^j k^kj, 



(14) 



describing two cones, which is a general situation for lev- 
els crossing [23|. However, in the case considered, we can 
be more specific about these cones. Because any vector 
in the k x , k y plane compatible with the symmetry C% v is 
identically equal to zero, and any tensor of rank two com- 
patible with the symmetry is proportional to the unity 
tensor, the dispersion law is just 



£~±|k| 



(15) 



and the cones are circular, with the axis perpendicular 
to the k x , k y plane. 

The important role played by discrete symmetries in 
protecting a fc-linear dispersion in graphene was pointed 



out in Ref. [24J. Ref. [25j was an important precedent 



in applying group theory methods to graphene. The ap- 
pearance of massless Dirac fermions under conditions of 
hexagonal symmetry was considered in Ref. [3]. Group 
theory was used to derive an invariant expansion of the 
Hamiltonian for electron states near the K points of the 
graphene Brillouin zone in Ref. [27l |. The influence 
of stress on the bands merging was analyzed in Refs. 
HI, HI- The influence of spin-orbit interactions on the 
band structure of graphene was studied in Ref. [3fl | . 
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